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ABSTRACT 

We probe the statistical nature of the primordial density field by measuring correla- 
tions between the power in adjacent Fourier modes. For certain types of non-Gaussian 
field, these fc-space correlations would be expected to be more extended than for a 
Gaussian field, providing a useful discriminatory test for Gaussian fields. We apply 
this test to the combined QDOT and 1.2-Jy IRAS galaxy survey and find the observed 
density field to be in good agreement with having Gaussian density fluctuations for 
modes with k < 0.1 h Mpc -1 . From this result we are able to set quantitative limits 
on a class of possible non-Gaussian distributions - the product of a Gaussian field 
with an independent stochastic field. The maximum sensitivity is to modulations of a 
Gaussian field with coherence scales of ~ 30 h^ 1 Mpc and the rms modulation on this 
scale cannot greatly exceed unity. We discuss the improvements to this limit likely to 
be set by future surveys. 



1 INTRODUCTION 

A number of theories now set out to explain the initiation 
of density perturbations in the early universe. Since most of 
these ideas invoke new physics, with the prediction of addi- 
tional fundamental fields, it is of interest to any cosmologist 
to narrow down the alternatives. One testable prediction 
that separates these theories is the statistics of the initial 
perturbations - some predict these to be non-Gaussian in 
form (most topological defect theories, and some versions of 
inflation), and others predict that the density field should 
be Gaussian (most types of inflation). 

Previous work in this area has concentrated mainly on 
three statistical approaches: (1) comparing the moments 
of the density probability distribution function found from 
galaxy redshift surveys with those predicted from a Gaussian 
distribution, (e.g. Saunders et al. 1991; Gatzanaga 1992; 
Nusser et al. 1995); (2) using topology to estimate the genus 
for the observed density-field contours (for instance Coles et 
al. 1993, Gott et al., Moore et al. 1992; Park et al. 1992; Vo- 
geley et al. 1994); (3) using N-body simulations to compare 
the evolution and clustering properties of galaxies, starting 
from different non-Gaussian and Gaussian initial conditions 
(e.g. Moscardini et al. 1991, Weinberg & Cole 1992.) 

Much of the work on higher order moments has con- 
centrated on measuring statistics of the density field in 
the quasi-linear regime, and comparing the observed non- 
Gaussian characteristics with what would be expected for 
the gravitational evolution of initially Gaussian fluctuations. 
In the above analyses the density field has been found to 
be consistent with a distribution that was initially Gaus- 
sian. However, in all this work it is necessary to smooth the 
galaxy distribution with some filter, averaging the effect of 



different regions of space. The central limit theorem then 
suggests that there is the danger that the appearance of 
Gaussian statistics will always be produced, whatever the 
underlying distribution. 

An alternative approach adopted by Feldman, Kaiser, 
and Peacock (1994; FKP) is to look at the density field in 
Fourier space, where modes are separated out on the basis of 
scale, making it easy to probe the linear regime. FKP pro- 
posed two possible tests for a Gaussian density field. Both 
make use of the property that at high resolution in k space, 
independent Fourier modes have power values that fluctuate 
about the mean power for that particular scale. The first test 
was to look at the one-point distribution of these power fluc- 
tuations, which is predicted to be exponentially distributed 
in the Gaussian case. FKP found that for the QDOT l-in-6 
survey, the power distribution was in good agreement with 
the Gaussian prediction. However, Fan & Bardeen (1995) 
have subsequently argued that such one-point pdf tests have 
little discriminatory power: for a large enough sample, most 
types of non- Gaussian field are also expected to have an ex- 
ponential distribution as a consequence of the central limit 
theorem. The second test suggested by FKP was to look 
at correlations between the power fluctuations which, for a 
Gaussian field, will be a function of the selection function 
only. This is less susceptible to the swamping effects of the 
central limit theorem, and it is the test that we investigate 
in this paper. 



2 POWER CORRELATIONS 

The general idea behind looking for correlations between 
power fluctuations is that for any infinite density field the 
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Fourier modes are independent i.e.: 
{S kl <5k 2 ) = (27r) 3 P(fci) * D (ki -ka), 



(1) 



where 5D(ki — k2) is the Dirac delta function and angle 
brackets here, and throughout the rest of this paper denote 
the ensemble average. For a finite sample of the density field, 
however, this mode independence is lost. This is because 
the finite region is effectively a product of the true infinite 
field, <5 mflnlte j an d a mask, m. In Fourier space the observed 
density modes (5° bsorvod are then a convolution of the 'true' 
density modes, with the Fourier transform of the mask. 



= 6? 



where <5k is given by 
<5k 



5{r) e <k r d 3 r, 



(2) 



(3) 



and S(r) = [p(r) — p(r)]/p(r) is the overdensity at position r. 
This has the effect of mixing information between the modes 
over a scale comparable to one over the fe-space width of 
the mask. Now, for a Gaussian field, the following relation 
applies (see FKP appendix B): 



|(«5 kl «5 k2 )| 2 = (P(kO P(k 2 ) ) - P(kO P(k 2 ) 
= (<5P(kO <5P(k 2 ) ) 



(4) 



where P(k), the ensemble-average power, comes from equa- 
tion (1); P(k) is the power in a single mode: 



P(k) = \6 k [ 



(5) 



and SP (k) = P (k) — P (k). This relation is dependent on 
two properties of Gaussian fields - a normal distribution 
for the one-point pdf, and the independence of the Fourier 
modes (i.e. the modes should have uncorrelated phases). For 
a finite sample of a Gaussian field, there will thus be corre- 
lations between the different power modes, with a two-point 
function for the power fluctuations which is the square of 
the two-point function for the Fourier amplitudes. For non- 
Gaussian fields in which intrinsic phase correlations exist, 
this proof does not apply, and we consider in Section 4 spe- 
cific non-Gaussian models in which the power correlations 
are broader than the Gaussian prediction. This does not 
prove that all conceivable non-Gaussian fields could be de- 
tected in this way, but it does show that the correlations 
of the fluctuating power field give a necessary condition for 
the field to be considered Gaussian. Figure 1 shows how the 
measured field appears for the IRAS-galaxy data described 
in Section 3. 

Quantitatively, we can define a power correlation func- 
tion, as follows: 

f (Ak) = ( [P(k) - P(fc)] [P(k - Ak) - P(fc - Afc)] ) 
( [P(k) - P(k) ] 2) 

This is simply |(5ki<5k 2 )| 2 as in equation (4), but normalised 
so that £p(0) = 1. We now show how to evaluate this func- 
tion for practical datasets. 

For a galaxy survey FKP defined the 'weighted galaxy 
fluctuation field', S(r) as 



S(r) 



TO ( r )[ n g( r ) ~ <*n a (r)] 
[J #rn 2 (r) w 2 (r) ] V 2 






(7) 
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Figure 1. Illustrating the fluctuating power field from the 
IRAS data described in Section 3. Equal area projections of 
power from shells of constant |k| are shown (Figure la at 
k = 0.1ft Mpc" 1 , and figure lb at k = 0.15 h Mpc" 1 ). Black 
corresponds to low, and white to high power. The panels have 
width Afc = 0.25 ft Mpc^ 1 . While the power is not constant for 
a particular value of |k|, there is a certain coherence scale over 
which there are significant correlations in power. As expected, this 
coherence scale does not appear to change with differing values 
of |k|, and it forms the basis for a test for Gaussianity. 



where n g (r) = E; 5n(r — rj) with r; being the galaxy posi- 
tion vectors. The smooth background density is subtracted 
via the analogous number density n 3 , which applies to a syn- 
thetic catalogue with number density 1/q times that of the 
real catalogue. n(r) is the expected mean density of galaxies 
given the angular and luminosity selection criteria, and w(r) 
is a weighting function designed to minimise the variance in 
power by favouring distant galaxies while the shot noise is 
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not a dominant contributor to the power. We have used the 
weighting function from FKP, which is: 



w(r) = [l + n(r)P(fc) 



(8) 



If this survey has an underlying Gaussian distribution, 
then FKP showed that the power correlation function would 
be of the form: 



Hp (Ak) = 



P(fc)Q(Ak) + S(Ak) 



(9) 



[P(k)+S(0)}* 

provided the mean power, P can be taken as constant over 
the width of Q. This expression contains two functions of 
the mean density: 



Q(k) 
and 

S(k) = 



/An 2 (r) w 2 (r) e ikr 
J d 3 r n 2 (r) w 2 (r) 

(1 + a) / d 3 rn(r) w 2 (r) e i 
J d 3 r n 2 (r) ui 2 (r) 



(10) 



(11) 



The idea now is to estimate £p via equation (6) and to com- 
pare this with the Gaussian prediction quoted in equation 
(9). 



3 APPLICATION TO IRAS GALAXIES 
3.1 Method 

We have found the power correlation function, £p for a com- 
bined dataset of ~ 4500 IRAS galaxies, consisting of the 
QDOT l-in-6 0.6-Jy redshift survey (see Efstathiou et al., 
1990; Lawrence et al. 1996) and the Berkeley 1.2-Jy red- 
shift survey (Fisher et al. 1993, Fisher et al. 1995). The 
power in each Fourier mode, (= <5k| 2 ) was found following 
the procedure described in FKP (see also Tadros & Efs- 
tathiou 1995). We have used a direct Fourier transform of 
the survey, choosing ~ 6000 random orientations in the k- 
space shell over which to perform the integration. This was 
to sample power modes close enough together in fc space to 
determine the shape of the correlation function accurately. 
For an estimator of P(k), we have averaged P(k) over these 
different orientations in the fc-space shell: 



P(k) 



M N 



(12) 



=1 i=i 



£p was then evaluated using the assumption that the en- 
semble average can be replaced by an average over different 
orientations in k space. This assumption holds for a Gaus- 
sian field, which is the hypothesis we will be testing against. 
There is a small amount of anisotropy introduced by the 
angular part of the mask, so in comparing £p(Afc) with its 
theoretical prediction, we have also averaged Q, and S over 
the different orientations in the fc-space shell. 

The normalisation of £p depends on the correct evalu- 
ation of the mean power, P, which can lead to problems if 
there are not enough independent patches in fc space to get 
an accurate estimate. This can be overcome by evaluating 
£p on a shell of fixed |k|, for which the mean power is a 
constant. Then by requiring that £p — > for large Afc, we 
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Figure 2. The power correlation function for the QDOT + 1.2-Jy 
survey. The solid line is the theoretical prediction for this survey 
if the underlying distribution is Gaussian. Error bars are based 
on Monte Carlo simulations of the deviations from the theoretical 
form for Gaussian surveys of the same size as the QDOT + 1.2-Jy. 

can find the actual value for P, and renormalise accordingly. 

The values of |k| used were chosen experimentally using 
the criteria that the scales probed should be in the linear 
regime, but not so large that shot noise dominates the signal. 
The fc-space shell also needed to have large enough radius 
for there to be enough independent coherent patches over 
which to find an average power correlation function. 

The galaxy coordinates that we use are in redshift 
space, and so peculiar velocities may affect the predicted 
shape of £p. Empirically we find this to be a small effect, 
in the sense that £p appears to be the same the radial and 
transverse directions, and we will examine it in detail else- 
where. 

3.2 Results 

Figure 2 is a plot of the power correlation function for the 
QDOT and 1.2-Jy survey with the prediction for an un- 
derlying Gaussian distribution. The error bars are based on 
expected deviations from the curve for a catalogue that has 
an underlying Gaussian distribution. These have been de- 
termined by making real space Monte Carlo realisations of 
the galaxy catalogue, and finding the spread in the shape of 
the correlation functions. The realisations are generated by 
Poisson sampling from a given underlying probability dis- 
tribution, giving all catalogues the same selection function, 
angular mask, and power spectrum as the real catalogue, 
with a Gaussian distribution on large scales. The density 
fluctuations are generated as lognormal realisations (Coles 
& Jones 1991) to give non-linear features on smaller scales. 

The figure shows that the QDOT+1.2 -Jy survey ap- 
pears to lie within the error bars predicted for an initial 
Gaussian distribution. Quantitatively, the goodness of fit 
can be characterised by performing a x 2 -type test on the 
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Figure 3. Predicted shapes of the correlation function for a prod- 
uct of Gaussian field: <5(r) = 5c(r)[l + r?(r) ] for varying forms of 
the power spectrum of rj, {\r]k\ 2 )- in the upper graph, the mod- 
ulating scale, R c is fixed at 50/i _1 Mpc, and the different lines 
represent different values of a 2 : From the bottom up, the lines 
are for (1) zero amplitude of modulation (i.e. a Gaussian field); 
(2) a 2 = 0.1; (3) a 2 = 1.0; and (4) (the uppermost) a 2 = 10.0 
The lower graph has a fixed amplitude, a 2 of 3.0, with varying 
modulating scales. From the bottom up, line (1) is for the Gaus- 
sian prediction as in the graph above; (2) R c = 150 h^ 1 Mpc; (3) 
the dashed line is for R c = 30 h' 1 Mpc; (4) R c = 100 h" 1 Mpc; 
and (5) (the uppermost line) for R c = 50 hr 1 Mpc. 

data: 



(13) 



where Oi is the set of observed data points, d is the equiva- 
lent set expected for a Gaussian distribution, and Oi is the 
standard deviation for the observed data points, which we 
have taken to be the same as that found for the Gaussian 
case from the Monte Carlo simulations, ip 2 does not follow 
the conventional \ distribution, since the data points them- 
selves are correlated, making the actual number of degrees 
of freedom fewer than the number of data points. 

The probability distribution for ip 2 has been determined 
from a frequency plot of ip 2 values found from the 100 Monte 
Carlo simulated galaxy catalogues mentioned above. We find 



a probability, p(ip 2 > V'moasurod) of 63%, (for the null hypoth- 
esis that the data set follows the Gaussian prediction). This 
is in good agreement with the Gaussian prediction, but how 
does it compare with what any non-Gaussian model might 
predict? 



4 LIMITS ON NON-GAUSSIAN MODELS 

Since the power correlation function for the Gaussian case is 
a function of the selection function only, one particular class 
of non-Gaussian models can be analysed with relative ease. 
This class is the product of a Gaussian field with another 
stochastic field, uncorrelated with the former: 



6(r) =$ G (r)[l + »j(r)] 



(14) 



This has been suggested as a possible alternative form for the 
density field, both on empirical grounds by Peebles (1983), 
and in some inflationary scenarios involving multiple scalar 
fields (e.g. Yi & Vishniac 1993). Locally, this field looks 
Gaussian, but on larger scales it is modulated so that the 
amplitude of fluctuations varies in different parts of space 
giving rise to quiet and noisy regions. This is equivalent to 
a field with zero skewness, but non-zero kurtosis. The one- 
point pdf is a simple test for this form of non-Gaussian be- 
haviour; however a simple transformation of the field would 
restore a Gaussian looking pdf, but it would still be modu- 
lated in a way which our method could detect. 

We can predict the form of £p by treating the (1 + rf) 
part of the field as part of the selection function, so that 



nw — > nw(l + ij) 

in equations (10), and (11). Q then becomes: 
Q' = Q + 2 Q * 77k + Q * r? k * r? k , 
and S: 

S' = S + S * r? k 



(15) 
(16) 
(17) 



We can now use equation (9) to obtain the modified corre- 
lation function: 



& (Ak) 



P{k)Q' (Ak) + S' (Ak) 
[P(k)Q>(0)+S> (0)] 2 



(18) 



Substituting in the terms for Q' and S' from equations (16) 
and (17), and setting PQ' (0) + S' (0) = P' for ease of nota- 
tion, this becomes: 



P' 2 Cp (Ak) = (|PQ + S-| 2 ) 
+ (| (2PQ + 5)*r, k + PQ*r ?Ak *7 7Ak | 2 ) 
+ {(PQ + S)[ (2PQ + S)* r> Ak + PQ * r, A k * VAkT 
+ conj.) 



(19) 



Now, by noting that odd powers in r/u average to zero, and 
that (| / * f? k | 2 ) = (| / | 2 ) * (| r/ k | 2 ), where / is a function 
uncorrelated with rj, we can rewrite equation (19) as: 

P' 2 £ P (Ak) = i G P + (\ 2P(k) Q + S \ 2 )*( | r /Ak | 2 ) 
+ ( | P(k) Q * rjAk * rjAk | 2 ) (20) 
+ ( (P(k) Q + S) [P(k) Q * r/ Ak * r; Ak ]* + conj ) 
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and this can be written as: 

P' 2 ^(Ak) =$(Ak) 
+ ( \2P(k)Q (Ak) + S (Ak) | 2 ) * <|77Ak|> 2 
+ P(fc) 2 (!Q(Ak)|) 2 *(| J7Ak | 2 *|r /Ak | 2 ) 
+ 2( Re{P{k)Q (Ak) [P(k)Q (Ak) + S (Ak)]*}) x 



(21) 



/ 



(M 2 ) dk, 



where £p is the Gaussian contribution to the power corre- 
lation function. The effect of r] is to broaden the correlation 
function by an amount dependent on the power spectrum 
of the modulating field, thus acting like an unaccounted for 
part of the selection function. Since this non-Gaussian field 
is a product of a Gaussian with a modulating field, and the 
modulating part can be treated as part of the selection func- 
tion, we are still justified in our use of the approximation in 
equation (12). 

Figure 3 shows different shapes of the power correla- 
tion function, for a range of power spectra for r). We have 
parameterised this power spectrum by supposing that the 
stochastic process, ij, is white noise with an amplitude a 2 , 
and a Gaussian cut-off on scales R c , so that 



Pv = (M 2 )= (2 ^i? c )Vexp[-fc 2 .R c 2 ] 



and 



P„ d 3 k = a 2 



(22) 



(23) 



(27T) 

For the simple case of a small amplitude modulation in the 
absence of shot noise, equation (21) reduces to: 

{\Q'(Ak)\ 2 ) = {\Q + 2Q* Vk \ 2 } 
= exp[-Afc 2 /4a 2 ] 



+ 



"S2o?R 3 c a 2 
(1 + 4a 2 R 2 ) :i / 2 



(24) 



cxp[-Afc 2 i? 2 /(4a 2 i? 2 + 1)] 



where a is the effective fc-space Gaussian size of the survey. 
This expression is illustrated in figure 3. 

We are now in a position to constrain the amplitude 
and scale of modulation that could possibly be present in the 
observed density field. We have the measured value of ip 2 for 
the QDOT+1.2Jy survey from the Gaussian prediction, and 
can specify the largest allowed value of ip 2 at a given level of 
confidence. For a given deviation, we can find a set of values 
of a and R c for the power spectrum of r) which would give 
this level of departure from the Gaussian prediction. These 
contours are shown in figure 4. Present data allow a ~ 3 
modulation on ~ 30 h -1 Mpc scales at the 95% level. 

It is clear from figure 4 that the sensitivity to modula- 
tions is peaked around scales of ~ 30 ft -1 Mpc. From figure 
3, we see that this scale is set by the width of £p for the un- 
modulated field. Both larger and smaller scales give a smaller 
perturbation to £p; in the former case because there are too 
few independent samples within the sample volume. The lat- 
ter case is more subtle: although modulations characterised 
by a small length scale give rise to broad power correlations, 
the amplitude of the modulation is small compared with the 
amplitude of the Gaussian field. So fluctuations in the power 
modes persist for small Ak - as with a Gaussian field, and 
the effect of the modulation is only seen at large Ak. We are 
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Figure 4. Contour plot of ip 2 deviations in R c — a 2 param- 
eter space. The contours represent different ip 2 values for the 
deviation of £p from the Gaussian theoretical prediction. The 
outermost contour is the measured deviation, ip 2 of §p for the 
QDOT+1.2— Jy survey from the Gaussian form. This corresponds 
to a probability that ip 2 > 1 / ,2 ncasurod °f 63% if the survey is from 
a Gaussian distribution. The innermost contour indicates the re- 
gion of parameter space for (|%| 2 } that can be ruled out with 95% 
confidence as deduced from the form of £p for the QDOT+1.2-Jy, 
and the corresponding error bars. 



least sensitive to large Afe correlations because of the need to 
renormalise the correlation function due to the uncertainty 
in the mean power, P (see Section 3.1). 



5 SUMMARY AND OUTLOOK 

This preliminary exploration suggests that the use of power 
correlations can be a useful tool in allowing quantitative 
limits to be placed on certain types of non-Gaussian fields. 
Future work will include examining in more detail the effects 
of redshift distortions, and extrapolating the technique for 
use in larger redshift surveys. Since the error bounds scale 
roughly as volume -1 / 2 , we can expect to tighten the limits 
on a by about one order of magnitude with the next genera- 
tion of redshift surveys. We should also be able increase the 
sensitivity of this test to modulations on scales approaching 
600 h~ x Mpc. It will be interesting to compare the power of 
this test with alternative methods, such as the use of topol- 
ogy to detect non-Gaussian signatures. 
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